Ultracold Plasma Expansion in a Magnetic Field 
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We measure the expansion of an ultracold plasma across the field lines of a uniform magnetic 
field. We image the ion distribution by extracting the ions with a high voltage pulse onto a position- 
sensitive detector. Early in the lifetime of the plasma (< 20 /xs), the size of the image is dominated 
by the time-of-fiight Coulomb explosion of the dense ion cloud. For later times, we measure the 
2-D Gaussian width of the ion image, obtaining the transverse expansion velocity as a function of 
magnetic field (up to 70 G). We observe that the expansion velocity scales as B~^", explained by 
a nonlinear ambipolar diffusion model with anisotropic diffusion in two different directions. 

PACS numbers: 52.25.Xz, 52.55.Dy 



Plasma expansion in a uniform magnetic field is of in- 
terest in astrophysical, ionospheric, and laser-produced 
plasma applications. The presence of a magnetic field 
during expansion can initiate various phenomena, such 
as plasma confinement and plasma instabilities [ij. Ul- 
tracold plasmas (UCPs), formed by photoionizing laser- 
cooled atoms near the ionization limit, have system pa- 
rameters many orders of magnitude away from tradi- 
tional laser-produced plasmas, with electron and ion tem- 
peratures on the order of meV, or even /leV, and densities 
of 10^ to 10^" cm^"^. UCPs thus provide a testing ground 
to study basic plasma theory in a clean and simple sys- 
tem, and we need fields of only tens of Gauss to observe 
significant effects on the expansion dynamics (for our sys- 
tem the electrons are magnetized, while the ions are not). 
All previous studies of ultracold plasma expansion, both 
experimental [2, y, Ij, l5| and theoretical \§., J., 3., M i focus 
on free expansion in the absence of magnetic fields. 

In this work, we present the first measurement of UCP 
expansion in a magnetic field. By extracting the ions 
with a high voltage pulse onto a position-sensitive detec- 
tor, we image the ion distribution of an UCP expand- 
ing in a uniform magnetic field. Early in the lifetime 
of the plasma, the image size is dominated by the time- 
of-fiight Coulomb explosion of the dense ion cloud. At 
about 20 /IS the image size is at a minimum and then 
afterwards increases, reflecting the true size of the ex- 
panding plasma. The expansion is self-similar, as the ion 
cloud maintains a Gaussian density profile throughout 
the lifetime of the plasma. By 2-D Gaussian fitting of 
the ion image, we obtain the width transverse to the ap- 
plied magnetic field. In the absence of a magnetic field, 
the plasma expansion velocity at different initial electron 
temperatures (Te) matches the result obtained by mea- 
suring the plasma oscillation frequency [2, llO|. As we 
increase the field up to 70 G, we find that the transverse 
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FIG. 1: Experimental setup and two-photon excitation pro- 
cess for an UCP in a magnetic field, (a) two-photon excitation 
process: one photon (red solid line, vr polarized) at 882 nm 
drives the 6s[3/2]2 (m = 0) — > 6p[5/2]3 (m = 0) transition, 
and the other (10- ns pulse) at 514 nm ionizes the atoms in the 
6p[5/2]3 (m — 0) state. The dashed lines are the repumping 
beams (both blue and red of the 6s[3/2]2 ^ 6p[5/2]3 tran- 
sition, also TT polarized); (b) experimental setup for imaging 
the ions onto the MCP/phosphor screen. 



expansion velocity decreases, roughly scaling as B~^/^. 
This field dependence is well explained by an ambipolar 
diffusion model which involves anisotropic diffusion in 
two different directions: the diffusion rate is almost un- 
affected in the direction along the magnetic field, while 
it is reduced in the direction normal to the field. A crit- 
ical feature of this diffusion model is its dependence on 
density and temperature which, for an expanding UCP, 
are time dependent. 

Our production of UCPs proceeds as detailed in 
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Metastable Xenon atoms are cooled and trapped in a 
magneto-optical trap (MOT). The neutral atom cloud 
has a temperature of about 20 /xK, total number of about 
2 X 10^, peak density of about 2 x 10^° cm"'^, and a 



Gaussian spatial density distribution with an rms radius 
of about 280 /jLir. The plasma is then produced by a two- 
photon excitation process (figure la), ionizing up to 30% 
of the atoms. One photon for this process is from the 
cooling laser at 882 nm, and the other is from a pulsed 
dye laser at 514 nm (10-ns pulse). We control the ioniza- 
tion fraction with the 514-nm intensity, while the initial 
electron energy is controlled by tuning the 514-nm pho- 
ton energy with respect to the ionization limit. In the 
absence of a magnetic field, we can simply use the MOT 
beams to drive the 6s [3/2] 2 — > 6p[5/2]3 transition; with 
the addition of a magnetic field, the 6s[3/2]2 and 6p[5/2]3 
states split into 5 and 7 Zeeman sub-levels respectively, 
shifting the transitions out of resonance. To maintain a 
reasonable excitation fraction, we need to optically pump 
the atoms into the m = magnetic sublevel of the ground 
state and then drive the 6s [3/2] 2 (m = 0) — > 6p[5/2]3 (m 
= 0) transition, which will be unaffected by the magnetic 
field to first order. (We cannot turn on magnetic fields in 
a time short compared to the time that ions acquire ve- 
locity in the plasma due to electron pressure (^ 1 fis) and 
thus must photoionize in the presence of the field). It is 
difficult to accumulate a large m = population via opti- 
cal pumping for a J ^ J + 1 transition, so we developed 
a multiple frequency scheme that uses Zeeman shifts to 
address the m = and m = ± 1 sublevels (figure la). We 
use frequencies detuned -5, 0, 5 MHz respectively relative 
to the 6s [3/2] 2 (m = 0) ^ 6p[5/2]3 (m = 0) transition, 
and generate three 7r-polarized optical pumping beams. 
With this configuration, we can achieve UCP densities of 
~ 50% of the zero field value, approximately independent 
of magnetic field up to 100 G. 

The ionized cloud rapidly loses a few percent of the 
electrons, resulting in a slightly attractive potential 
for the remaining electrons, and it quickly reaches a 
quasineutral plasma state. It then expands with an 
asymptotic velocity typically in the 50-100 m/s range 
caused by the outward electron pressure [2|. By ap- 
plying a high- voltage pulse to the top grid (about 1.5 
cm above the plasma) and with accelerating voltages on 
the middle and front grids (figure lb), we image (mag- 
nification of 1.3) the ion distribution of the UCP on a 
position-sensitive microchannel plate detector with phos- 
phor screen. The phosphor images are recorded with a 
CGD camera. The high-voltage pulse has an amplitude 
of 340 V, a width of 4.5 /zs, and a rise time of 60 ns. We 
use a modified square pulse generator [l2], which uses 
power FETs to fast switch a high voltage source. The 
magnetic field is generated by a pair of coils which are 
located outside the vacuum chamber, and is fully turned 
on before the creation of the UCP. Figure 2 is a typical 
result with an average of 8 images to increase the signal- 
to-noise ratio. Figure 2a is a false color plot of an ion 
image after 20 /xs of plasma evolution for Te{Q) = 100 K 
and B = 20 G, which fits well to a 2-D Gaussian profile 
(figure 2b). The ion images maintain a Gaussian profile 
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FIG. 2: (a) a false color 2D ion spatial distribution integrated 
over the third dimension of the UCP for t = 20 ^s, Te = 100 
K and B = 20 G; (b) the 2D Gaussian fittings (blue curves) 
of the ion image (a) along the x and y axis in the horizontal 
plane (red curves) ;(c) and (e) are the contour plots of the ion 
images at different delay times for G and 50 G respectively; 
(d) and (f) are the corresponding 2-D Gaussian fit of (c) and 
(e). All the units in (a)-(f) are in pixel number, and 1 pixel 
unit is about 150 /im. 



during the whole lifetime of the UCP (about 200 ^s) as 
shown in figure 2 c-f. By 2-D Gaussian fitting of the ion 
image, we extract the size transverse to the magnetic field 
at specific delay times. In figure 3, we plot the transverse 
sizes of UCPs in a field as a function of time. The mea- 
sured size of the image for early times is dominated by 
the Coulomb explosion of the ions, and reflects the overall 
Coulomb energy rather than the plasma size 10, ll3| . As 
the plasma size increases, the Coulomb explosion effect 
diminishes (note decreasing images at early times). After 
about 20 /iS, the measured image reflects the plasma size, 
and the images increase with increasing plasma time, as 
expected. As we increase the field, the transverse size 
increases more slowly with time, i.e., the slopes of the 
curves decrease, which indicates a slower expansion and 
magnetic confinement in the transverse direction of the 
UCP. For G at long times, the size does not linearly 
increase. This is in part because the size of the UCP is 
large enough to be affected by the four posts that secure 
our grids above and below the plasma. We also notice 
that the ion images have a flat top and even dip at late 
times of about 150-200 /is for the G case, the cause of 
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FIG. 3: Transverse sizes of UCPs in a magnetic field. The 
black curve (diamonds) is for G, the red curve (triangles) is 
for 20 G, and the blue curve (dots) is for 50 G. Early in the 
lifetime of the plasma, the size of the image is dominated by 
the Coulomb explosion of the dense ion cloud. 



By linear fitting to the transverse sizes vs. time after 
about 20 /iS (for G, only fitting the restricted linear 
region), we extract the transverse asymptotic expansion 
velocity of the UCP. In figure 4, we plot the transverse 
asymptotic expansion velocity as a function of magnetic 
field. The blue solid curve with square points is the ex- 
perimental result, the black dashed curve is a power-law 
fit of the data, and the red curve with triangle points 
is a simulation by assuming that the transverse expan- 
sion is due to ambipolar diffusion (discussed below). In 
the absence of a magnetic field and Te(0) = 100 K, the 
UCP expands with an asymptotic velocity of about 70 
ni/s caused by the outward electron pressure. This is the 
ion acoustic velocity for re(0). As the field increases, the 
transverse expansion velocity decreases, roughly scaling 
as B~^/^. At 50 G, it is about 17 m/s, about a factor of 
four smaller than that at G. 

Using the same technique, we extract the transverse 
asymptotic expansion velocity for different Te(0)s, shown 
in fig. 5. The red solid line with square points is the ex- 
perimental result for G, which agrees with the results 
obtained by measuring the plasma oscillation frequency 
(the black solid curve with circle points) ^. This vali- 
dates our imaging method for measuring the transverse 
expansion velocity and previous thechnique. The dashed 
lines are linear fits of the data above 60 K. For 50 G, 
the transverse expansion velocity as a function of Te(0) 
is similar to that at G except for a factor 3-4 smaller. 
The expansion velocity goes as Tg ' for initial electron 
temperature higher than 60 K, that is, the slopes of the 
dashed lines in figure 5 are about 1/2. 

The expansion of the UCP at the ion acoustic velocity 
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FIG. 4: The transverse asymptotic expansion velocity vs. 
applied magnetic field. The blue curve (squares) is the exper- 
imental result, the black dashed curve is a power-law fit, and 
the red curve (triangles) is a simulation using an ambipolar 
diffusion model that involves anisotropic diffusion in two di- 
rections. The error bars represent the la standard uncertainty 
resulting from the linear fits to the curves in figure 3. 
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FIG. 5: The transverse expansion velocity as a function of ini- 
tial electron temperature T'e(O). The red solid curve (squares) 
is the experimental result for G, which matches results ob- 
tained by measuring the plasma oscillation frequency (the 
black solid curve (circles)) [3. The brown solid curve (tri- 
angles) is for 25 G, and the blue solid curve (stars) is for 
50 G. The dashed lines are power-law fits above 60 K with 
powers near 1/2 (0.44, 0.54 and 0.63 from top to bottom, 
respectively) . 



can be understood in terms of electron pressure in a sim- 
ple hydrodynamic model [2], or alternatively as ambipo- 
lar diffusion, where oppositely charged particles (elec- 
trons and ions) diffuse together at the lower ion velocity. 
In a magnetic field, ambipolar diffusion allows charged 
particles to cross field lines, and produces a diffusion 
constant that scales as B^^, and a transverse decay size 
(expansion velocity) scaling of B~^ for conventional low 
temperature plasmas such as arc plasmas and strongly 



ionized plasmas in steady state (Q machine) [14|, in con- 
trast to our observed B~^'^ scaling. Another possibility 
is anomalous diffusion (Bohm diffusion) with a diffusion 
constant scaling as i?~^, which is intimately connected to 
plasma instabilities [lj,ll5|. However, observing the elec- 
tron signals as well as the single shot ion images with the 
magnetic fields, we do not see any signature of plasma 
instabilities. Our UCPs are rather novel plasma sys- 
tems, not only because of their low temperatures, but 
also because they are unconfined, inhomogeneous, and 
freely expanding into vacuum (the same properties hold 
for the UCP expansion along the magnetic field). There 
is also competition between adiabatic cooling and three- 
body recombination-induced heating leading to a time- 
dependent falling Te [18|. Given the evolution of the 
plasma, we cannot simply apply steady-state ambipolar 
diffusion theory, but must include the dynamics. In a 
strong magnetic field, the diffusion rate of a plasma be- 
comes anisotropic; the diffusion rate is almost unaffected 
in the direction along the magnetic field, while it is re- 
duced in the direction normal to the field. For our UCP 
expansion in the magnetic field, we assume ambipolar 
diffusion in both directions (along the field and normal 
to the field) and write the diffusion equation as follows by 
separating the diffusion rate into two terms in cylindrical 
coordinates: 
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where Di^ is the transverse ambipolar diffusion rate with 
-Dj_ = pi^ei oc „2"°i/2 ; Pe IS the clcctron gyroradius and 
Vei is the electron-ion collision rate, tii is the character- 
istic time with T|| = a{t)/vs{t), a(t) and Vs{t) are the 
plasma size and the expansion velocity along the mag- 
netic field respectively. The first term in the right hand 
side of eq.l corresponds to the diffusion normal to the 
field, and the second term is the diffusion along the mag- 
netic field. Normalizing all the variables in cq. 1 to their 
initial values, we get: 
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where h = -^,t = 4-, 

"0 ' to ' 
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with io = ■1°'° and run = -^. 

For simplicity, we assume no spatial temperature gra- 
dient for our UCP, and that the density profiles in both 
directions remain Gaussian throughout the plasma life- 
time. In order to solve eq. 2, we also need to know Te{t) 
in the magnetic field. Although we might expect changes 
in the Tg evolution in the presence of a large magnetic 
field compared with free expansion, we did not observe 
any significant difference in the three-body recombina- 
tion rates, which can be used to extract Te{t) [l8|. We 
therefore assume that Te{t) is independent of magnetic 



field for the sake of our diffusion model (eq. 2). Elec- 
tron temperature evolution in an expanding UCP is a 
complicated problem and has been studied using various 
methods [5|, [la, [iTI, ll8| . Using Dopplcr broadening of the 
ion optical absorption spectrum together with numeri- 
cal simulation, Te (i) was observed to follow a self-similar 
analytic solution during the first 20 /is for Te(0) greater 
than 45K, that is Te = Te{0)/[l + t^ /rj^) ^. We obtain 
TM.) up to 70 /xs by using three-body recombination rate 
[l|, which gives Te ex t"^"^ for Teo = 3A'. We use the 
following empirical equation for Te{t) for our diffusion 
model: 



Te = Te(0)/(l+iV^||o)>i<20Ms 
Te = Te(20/is)(f/20)~^■^t > 20Ais 



(3) 



By numerically solving eq. 2 with T(,{t) given by cq. 
3, we find the transverse size evolution for different mag- 
netic fields from 5 G to 70 G. We extract the transverse 
expansion velocity as a function of magnetic field, shown 
as the red curve (triangles) in figure 4. The simulation re- 
sult is in good agreement with the B~^'^ scaling observed 
in the experimental data. Therefore ambipolar diffusion 
can explain the B-dependence when the expansion dy- 
namics are included. Using self-similiar arguments for 
the diffusion [l^, it can be shown that r oc i?"i/2t(i+'^)/4 
where Te oc t~ . By taking the derivative of r, we find an 
expansion velocity scaling as 5"^/^, as observed. Even 
if Te(t) in the presence of magnetic field is different from 
eq. 3 (different A), it will only alter the time-dependence 
of the transverse size, but not the B~^/^ scaling of the 
above self-similiar argument. This scaling of the trans- 
verse size (expansion velocity) is the result of the elec- 
tron density dependence of the tranverse ambipolar dif- 
fusion rate, while that of other low temperature plas- 
mas is density independent and yields the usual B~^ de- 
pendence of the characteristic transverse decay distance. 
Bohm diffusion also has the same B~^l^ dependence of 
the transverse size for our UCP systems if we assume 
D ^ oc T^jeB. However, the plasma will quickly stop ex- 
panding in a very short time (< 20 yus) if we assume only 
Bohm diffusion is involved, because the Bohm diffusion 
rate is proportional to Tg. The ambipolar diffusion rate 

1 /9 

is inversely proportional to Tg ' , and should dominate 
Bohm diffusion at later times which is consistent with no 
signature of plasma instabilities. 

In conclusion, we have developed a method to study 
the transverse expansion of an ultracold plasma in a mag- 
netic field by using a high voltage pulse to image the ion 
distribution onto a position sensitive detector. The ex- 
pansion velocity roughly scales as B~^'^ for magnetic 
field up to 70 G, and is in good agreement with a non- 
linear ambipolar diffusion model. 
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